Bounds on the maximum multiplicity 
of some common geometric graphs 



Adrian Dumitrescu* Andre Schulz^ Adam Sheffer* Csaba D. T6th§ 

September 27, 2011 
i—i ■ 

O 

<N 

Abstract 

We obtain new lower and upper bounds for the maximum multiplicity of some weighted and, 
\Q ' respectively, non-weighted common geometric graphs drawn on n points in the plane in general 

position (with no three points collinear): perfect matchings, spanning trees, spanning cycles 
(tours), and triangulations. 

(i) We present a new lower bound construction for the maximum number of triangulations a 
set of n points in general position can have. In particular, we show that a generalized double chain 
formed by two almost convex chains admits f2(8.65") different triangulations. This improves the 
bound Q(8A8 n ) achieved by the previous best construction, the double zig-zag chain studied by 

O . Aichholzer et al. 

(ii) We obtain a new lower bound of 0(12. 00 n ) for the number of non-crossing spanning 
, trees of the double chain composed of two convex chains. The previous bound, 17(10.42"), stood 

^ ' unchanged for more than 10 years. 

(iii) Using a recent upper bound of 30" for the number of triangulations, due to Sharir 
and Sheffer, we show that n points in the plane in general position admit at most 0(68.62") 
non-crossing spanning cycles. 

(iv) We derive lower bounds for the number of maximum and minimum weighted geometric 
graphs (matchings, spanning trees, and tours). We show that the number of shortest tours can 
be exponential in n for points in general position. These tours are automatically non-crossing. 
Likewise, we show that the number of longest non-crossing tours can be exponential in n. It was 
known that the number of shortest non-crossing perfect matchings can be exponential in n, and 
here we show that the number of longest non-crossing perfect matchings can be also exponential 
in n. It was known that the number of longest non-crossing spanning trees of a point set can 
be exponentially large, and here we show that this can be also realized with points in convex 
position. For points in convex position we obtain tight bounds for the number of longest and 
shortest tours. We also give a combinatorial characterization of longest tours, which yields an 
0(n log n) time algorithm for computing them. 



in 



x 



Keywords: Geometric graph, non-crossing property, Hamiltonian cycle, perfect matching, triangu- 
lation, spanning tree. 

'Department of Computer Science, University of Wisconsin-Milwaukee, USA. Supported in part by NSF grants 
CCF-0444188 and DMS-1001667. Email: dumitres@uwm.edu 

^Institut fur Mathematische Logik und Grundlagenforschung, Universitat Miinster, Germany. Partially supported 
by the German Research Foundation (DFG) under grant SCHU 2458/1-1. Email: andre.schulz@uni-muenster.de 

^School of Computer Science, Tel Aviv University, Tel Aviv 69978, Israel. Supported by Grant 338/09 from the 
Israel Science Fund. Email: sheffera@tau.ac.il 

^Department of Mathematics and Statistics, University of Calgary, Canada. Supported in part by NSERC grant 
RGPIN 35586. Research by this author was conducted at Tufts University and at Ecole Polytechnique Federale de 
Lausanne. Email: cdtoth@ucalgary.ca 



1 



1 Introduction 



Let P be a set of re points in the plane in general position, i. e., no three points lie on a common 
line. A geometric graph G = (P, E) is a graph drawn in the plane so that the vertex set consists of 
the points in P and the edges are drawn as straight line segments between the corresponding points 
in P. All graphs we consider in this paper are geometric graphs. We call a graph non-crossing if 
its edges intersect only at common endpoints. 

Determining the maximum number of non-crossing geometric graphs on n points in the plane 
is a fundamental question in combinatorial geometry. We follow common conventions (see for 
instance [24] ) and denote by pg(-P) the number of non-crossing geometric graphs that can be em- 
bedded over the planar point set P, and by pg(re) = max|p| =n pg(P) the maximum number of 
non-crossing graphs an n-element point set can admit. Analogously, we introduce shorthand nota- 
tion for the maximum number of triangulations, perfect matchings, spanning trees, and spanning 
cycles (i. e., Hamiltonian cycles); see Table HJ For example, tr(re) = 0(/(n)) means that any n- 
element point set can have at most 0(f(n)) triangulations, and tr(n) = 0,(g(n)) means that there 
exists some n-element point set that admits tt(g(n)) triangulations. 



Abbr. 


Graph class 


Lower bound 


Upper bound 


Pg(n) 


graphs 


0(41.18") LH [15] 


0(207.84") [16, 26] 


cf(n) 


cycle-free graphs 


fi(12.26 n ) [new, Theorem [2] 


0(164.49") [16, 26^ 


pm(n) 


perfect matchings 


n*(3 n ) [TS] 


0(10.07") [24 


st(n) 


spanning trees 


$1(12. 00 n ) [new, Theorem [2] 


0(146.37") [16, [26] 


sc(n) 


spanning cycles 


0(4.64") ]15J 


0(68.62 n ) [new, Theorem [3] 


tr(n) 


triangulations 


ft(8.65 n ) [new, Theorem [Q 


0(30") [26 



Table 1: Classes of non-crossing geometric (straight line) graphs, current best upper and lower bounds. 



In the past 30 years numerous researchers have tried to estimate these quantities. In a pivotal 
result, Ajtai et al. [2] showed that pg(re) = 0(c n ) for an absolute, but very large constant c > 0. The 
constant c has been improved several times since then. The best bound today is c < 207.84, which 
follows from a combination of a result of Sharir and Sheffer [26] with a result of Hoffmann et al. [16] . 
Interestingly, this upper bound, as well as the currently best upper bounds for st(n), sc(n), and 
cf (re), are derived from an upper bound on the maximum number of triangulations, tr(re). This 
underlines the importance of the bound for tr(re) in this setting. For example, the best known 
upper bound for st(n) is the combination of tr(n) < 30" [26] with the ratio sc(n)/tr(re) = 
O* (4.879") [16]; see also previous work [22j [231 EH ES]. To our knowledge, the only upper bound 
derived via a different approach is that for the number of perfect matchings by Sharir and Welzl [24J , 
pm(re) = 0(10.07"). 

So far, we recalled various upper bounds on the maximum number of geometric graphs in cer- 
tain classes. In this paper we mostly conduct our offensive from the other direction, on improving 
the corresponding lower bounds. Lower bounds for unweighted non-crossing graph classes were 
obtained in [U EJ [15]. Garcia, Noy, and Tejel [15] were the first to recognize the power of the 
double chain configuration in establishing good lower bounds for the number of matchings, trian- 
gulations, spanning cycles and trees. It was widely believed for some time that the double chain 
gives asymptotically the highest number of triangulations, namely 0*(8"). This was until 2006, 
when Aichholzer et al. pQ showed that another configuration, the so-called double zig-zag chain, 
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admits @*(V7T) = 0(8.48 n ) triangulating. The double zig-zag chain consists of two flat copies 
of a zig-zag chain. A zig-zag chain is the simplest example of an almost convex polygon. Such 
polygons have been introduced and first studied by Hurtado and Noy [18J . In this paper we further 
exploit the power of almost convex polygons and establish a new lower bound tr(n) = 0(8. 65 n ). 
For matchings, spanning cycles, and plane graphs the double chain still holds the current record. 

Less studied are multiplicities of weighted geometric graphs. The weight of a geometric graph 
is the sum of its (Euclidean) edge lengths. This leads to the question: how many graphs of a 
certain type (e. g., matchings, spanning trees, or tours) with minimum or maximum weight can be 
realized on an n-element point set. The notation is analogous; see Table 2. Dumitrescu [9] showed 
that the longest and shortest matchings can have exponential multiplicity, 2^ n \ for a point set in 
general position. Furthermore, the longest and shortest spanning trees can also have multiplicity 
of 2^( n ). Both bounds count explicitly geometric graphs with crossings; however these minima are 
automatically non-crossing. The question for the maximum multiplicity for non-crossing geometric 
graphs remained open for most classes of geometric graphs. Since we lack any upper bounds that are 
better than those for the corresponding unweighted classes, the "Upper bound" column is missing 
from Table [2 



Abbr. 


Graph class 


Lower bound 


pm min (n) 


shortest perfect matchings 


0(2 n / 4 ) [9J 


P m max(™) 


longest perfect matchings 


0(2 n / 4 ) [new, Theorem [4] 


Stxnhi (71) 


shortest spanning trees 


0(2 n / 2 ) [9j 


Stmax(w) 


longest spanning trees 


fi(2 n ) [new, Theorem [5] 




shortest spanning cycles 


0(2 n / 3 ) [new, Theorem [7] 


SC max (/l) 


longest spanning cycles 


0(2 n / 3 ) [new, Theorem [6] 



Table 2: Classes of weighted non-crossing geometric graphs: exponential lower bounds. 
Our results. 

(I) A new lower bound, 0(8.65"), for the maximum number of triangulations a set of n points can 
have. We first re-derive the bound given by Aichholzer et al. [lj with a simpler analysis, which 
allows us to extend it to more complex point sets. Our estimate might be the best possible for 
the type of construction we consider. 

(II) A new lower bound, 0(12. 00 n ), for the maximum number of non-crossing spanning trees a 
set of n points can have. This is obtained by refining the analysis of the number of such trees 
on the "double chain" point configuration. The previous bound was 0(10. 42 n ). Our analysis 
of the construction improves also the lower bound for cycle-free non-crossing graphs due to 
Aichholzer et al. from 0(11. 62 n ) to 0(12. 23 n ). 

(III) A new upper bound, 0(68. 62 n ), for the number of non-crossing spanning cycles on n points in 
the plane. This improves the latest upper bound of 0(70. 22 n ) that follows from a combination of 
the results of Buchin et al. [6] and a recent upper bound tr(n) < 30 n by Sharir and Sheffer |26| . 

(IV) New bounds on the maximum multiplicity of various weighted geometric graphs on n points 
(weighted by Euclidean length). We show that the maximum number of longest non-crossing 
perfect matchings, spanning trees, spanning cycles, as well as shortest tours are all exponential 
in n. We also derive tight bounds, as well as a combinatorial characterization of longest tours 

The O*, O* , fi* notation is used to describe the asymptotic growth of functions ignoring polynomial factors. 
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(with crossings allowed) over n points in convex position. This yields an 0(n log n) algorithm 
to compute a longest tour for such sets. 

The main results of this paper were presented at the 28th International Symposium on Theoret- 
ical Aspects of Computer Science (STACS) in March 2011, in Dortmund, Germany. An extended 
abstract (without all the proofs) appeared in the corresponding conference proceedings [13]. Since 
then we were able to further refine the analysis in two of our main theorems. In particular, we 
improved the bounds in Theorems [2] and El 



1.1 Preliminaries 

Asymptotics of multinomial coefficients. Denote by H(q) = —qlogq — (1 — q) log(l — q) the 
binary entropy function, where log stands for the logarithm in base 2. (By convention, OlogO = 0.) 
For a constant < a < 1, the following estimate can be easily derived from Stirling's formula for 
the factorial: 

' H ) = ®{n- l l 2 2 H ^ n ), (1) 
an J 

We also need the following bound on the sum of binomial coefficients; see [5] for a proof and |1Q(, [TT] 
for an application. If < a < ^ is a constant, 

k<an 

' n 



£ ( ) <2 H ^ n . (2) 



k=0 

Define similarly the generalized entropy function of k parameters ax, . . . ,a k , satisfying 

k 



/ J aj = l, ai,... ,afe > 0, (3) 



8=1 



H k (ai, ...,a k ) = -y^ j a i log on 



i=l 



Clearly, H{q) = H2(q, 1 — q). Recall, the multinomial coefficient 



n \ n\ 



n 1 ,n 2 ,...,n k J n x !n 2 ! • • • n fe ! 



where Y2i=x n i = n i counts the number of distinct ways to permute a multiset of n elements, k of 
which are distinct, with n^, i = 1, . . . , k, being the multiplicities of each of the k distinct elements. 

Assuming that ni = a^n, i = 1, . . . , k, for constants a±, . . . , a k , satisfying (|3|), again by using 
Stirling's formula for the factorial, one gets an expression analogous to ([T|): 

(k \ n 

f[ ] = Q{ n -^l 2 ) ■ 2**(«i.-,«i> ( 4 ) 
i=l / 



Notations and conventions. For a polygonal chain P, let |P| denote the number of vertices. If 
1 < ci < C2 are two constants, we frequently write VL*(cV;) = fi(c"). A geometric graph G = (V,E) 
is called a (geometric) thrackle, if any two edges in E either cross or share a common endpoint; see 
e.g., [21]. 
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2 Lower bound on the maximum number of triangulations 



Hurtado and Noy [18] introduced the class P(n, k r ) of almost convex polygons. A simple polygon is 
in this class if it has n vertices, and it can be obtained from a convex r-gon by replacing each edge 
with a "flat" reflex chain having k interior vertices (thus there are r reflex chains altogether), such 
that any two vertices of the polygon can be connected by an internal diagonal unless they belong 
to the same reflex chain. For example, P(n, r ) is the class of convex polygons with n = r vertices. 
Note that, for r > 3 and k > 0, every polygon in P(n, k r ) has n = r{k + 1) vertices, r of which are 
convex. 

Following the notation from [18], we denote by t(n, k r ) the number of triangulations of a polygon 
P(n,k r ). According to [THJ Theorem 3], 

tin, n = e ((!+*/?) " ■ tin, o»>) = e ((* + *) r • 4*")) = ((<* + *)* . 2 )") . 

In particular, 

• for k = 1, t(n, V) = G((2 v / 3) n ) = G(v / 12 n ). This estimate was used by Aichholzer et al. pQ 
to show that the double zig-zag chain has Q(8A8 n ) triangulations. 

• for k = 2, t(n,2 r ) = 9((2 5 / 3 ) n ). 

• for k = 3, t{n, 3 r ) = 9((5 1 / 4 • 2) n ). 

• for k = 4, t(n,4 r ) = 9((6 1 / 5 • 2) n ). 

Analogously we define the class P + (n, k r ) of almost convex (polygonal) chains. An x-monotone 
polygonal chain in this class has n = r(k + 1) + 1 vertices, and it can be obtained from an x- 
monotone convex chain with r + 1 vertices by replacing each edge with a "flat" reflex chain with k 
interior vertices, such that any two vertices of the chain can be connected by a segment lying above 
the chain unless they are incident to the same reflex chain. See Fig. [1] for a small example. To 
further simplify notation, we denote by P + (n,k r ) any polygonal chain in this class; note they are 
all equivalent in the sense that they have the same visibility graph. Note that the simple polygon 
bounded by an almost convex chain P + (n,k r ) and the segment between its two extremal vertices 
has Q*(t(n,k r )) triangulations. 




Figure 1: Two (flat) mutually visible copies of P + (19, 2 6 ) with opposite orientations that form £>(38,2 12 ). 
Two consecutive hull vertices of P + (19, 2 6 ) with a reflex chain of two vertices in between are indicated in 
both the upper and the lower chain. 

In establishing our new lower bound for the maximum number of triangulations, we go through 
the following steps: We first describe the double zig-zag chain from [1J in our framework, and re- 
derive the Q*(\/72™) bound of pQ for the number of its triangulations. Our simpler analysis extends 
to some variants of the double zig-zag chain, and leads to a new lower bound on tr(n) = 17(8.65"). 
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Two x-monotone polygonal chains L and U are said to be mutually visible if every pair of points, 
p € L and q £ U, are visible from each other (that is, the segment pq crosses neither U nor L). 

Let us call D(n,k r ) the generalized double chain of n points formed by the set of vertices in 
two mutually visible x-monotone chains, each with n/2 = r(k + 1) + 1 vertices, where the upper 
chain is a P + (n, k r ) and the lower chain is a horizontally reflected copy of P + (n, k r ), as in Fig. [TJ 
Generalized double chains are a family of point configurations, containing, among others, the double 
chain and double zig-zag chain configurations. In particular, D(n, l r ) is the double zig-zag chain 
used by Aichholzer et al. pp. 



Theorem 1. The point set D(n,3 r 
sequently, tr(n) = $1(8.65") . 



with n = 8r + 2 points admits Sl(8.65 n ) triangulations. Con- 



Proof. The following estimate is used in all our triangulation bounds. Consider two mutually visible 
flat polygonal chains, L and U, with m vertices each (L is the lower chain and U is the upper chain). 
As in the proof of [151 Theorem 4.1], the region between the two chains consists of 2m — 2 triangles, 
such that exactly m — 1 triangles have an edge along L and the remaining m — 1 triangles have an 
edge adjacent to U. It follows that the number of distinct triangulations of this middle region is 



2m - 2 
m — 1 



Qim- 1 / 2 -4 r ' 



(5) 



The old 17(8.48"") lower bound in a new perspective. We estimate from below the number 
of triangulations of D(n,l r ) as follows. Recall that \L\ = \U\ = n/2 = 2r + 1. Include all edges 
of L and U in any of the triangulations we construct. Now construct different triangulations as 
follows. Independently select a subset of a\r short edges of conv(f7) and similarly, a subset of a.\r 
short edges of conv(L). Here a% £ (0, 1) is a constant to be chosen later. According to (UJ), this 
can be done in 

' r \ = Gfr- 1 / 2 • 2 H ^>) 
K axr) 

ways in each of the two chains. Include these edges in the triangulation. Observe that after adding 
these short edges the middle region between the (initial) chains L and U is sandwiched between 
two mutually visible shorter chains, say L' C L and U' C U, where 



\L'\ = \U'\ =2r- a x r 



Qi)r. 



(6) 



Triangulate this middle regions in all of the possible ways, as outlined in the paragraph above ([5]). 
Let denote the total number of triangulations of D(n, l r ) obtained in this way. By the above 
estimate, we have t(n, V) = 0((2\/3) n ). Combining this with © and ©, 



N = Q* 



n* 



[2\/3) 2r 2 H( - ai ^ r 4( 2 ~ Ql ) r 



2 2 . 3 . 2 (2-ai) 2 ^(«l) 



2r 



n* 



2<i— ai+H(ai) . 3 



n/2 



n* (a n ) 



where 



2 4— ai+H(oti) . 3 



(1/2) 



By setting ai = 1/3, as in [lj, this yields a = 6^/2 = 8.485. . ., and N = ft*(8.485 n ) = fi(8.48 n ). 

Applying a similar analysis for a generalized double chain with reflex chains of length 3 implies 
Theorem [TJ To simplify the presentation we describe the next step using almost convex polygonal 
chains with reflex chains of length 2. 
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The next step: using D(n, 2 r ). Notice that in this case n = 6r + 2. In the upper chain U, each 
reflex chain of two points together with the two hull vertices next to it form a convex 4-chain, as 
viewed from below. In each such 4-chain, independently, we proceed with one of the following three 
choices: (0) we leave it unchanged; (1) we add one edge, so that the chain length is reduced by 1 
(from 4 to 3); (2) we add two edges, so that the chain length is reduced by 2 (from 4 to 2). We refer 
to these changes and corresponding length reductions as reductions of type 0, 1 and 2, respectively. 
For i = 0, 1,2, let en count the number of distinct ways such reductions can be performed on a 
convex 4-chain of U. See Fig. [2jleft). Clearly, we have ao = 1 and a\ = a 2 = 2. 



Figure 2: Left: Length reductions for the chains of D(n, 2 r ): two of type 1 and two of type 2. Right: Length 
reductions for D(n, 3 r ): three of type 1, five of type 2, and five of type 3. 

For i = 0, 1, 2, make c^r reductions of type i in U, for some suitable (constant) values to be 
determined, with 



y^ai = l, a ,ai,Q 2 > 0, 



i=l 



According to (j3J), choosing one of the three reduction types for each of the 4-chains in U can 
be done in 



distinct ways. Once a reduction type i has been chosen for a specific 4-chain, it can be implemented 
in m distinct ways. It follows that the number of distinct resulting chains U' is 



0* / 2^3("o,"i I "2)r j JJ a Q i r = @* ( 2 H3 ^ a °' ai ' a2 ^ r ■ 2 aiT ■ 2 



So there are that many distinct chains U' which will occur. Proceed similarly for the lower 
chain L. Observe that once the ai are fixed, the two resulting sub-chains U' C U and L' C L have 
the same length 

\L'\ = \U'\ = 3r - air - 2a 2 r = (3 - a x - 2a 2 )r. (7) 

Triangulate the middle part (between U' and L') in any of the possible ways, according to ([5]). Let 
denote the total number of triangulations obtained in this way. Recall the estimate t(n,2 r ) = 
G((2 5 / 3 ) n ), which yields t(3r,2 r ) = 6((2 5 / 3 ) 3r ) = G(2 5r ). By © and ©, we obtain 



n = n* 
= 0* 



(2^) 2 • 2 2 ^ 3 ( a °' Ql,ct2 ) r • 2 2ctir • 2 2ct2r • 4-( 3 ~ ai ~ 



■2ct2)r 



where 



25 _ 2^3( a 0,c*i,cr2) . 2 ai . 2 Q2 • 2^ 3 ~ Ql ~ 2 ° 2 ) 
28— 02+^3 (ao, 01,02) 



2r 



n* (a n ) , 



1/3 
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The valued a = a x = 0.4 and a 2 = 0.2 yield a = 8.617. . ., and N = 0*(8.617 n ) = 0(8.61"). 

The new 0(8. 65 n ) lower bound: using D(n,3 r ). We only briefly outline the differences from 
the previous calculation. We have n = \D(n,3 r )\ = 8r + 2. For i = 0,1,2,3, let a, count the 
number of distinct ways reductions can be performed on a convex 5-chain of U. As illustrated in 
Fig. fright), we have ao = 1, ao = 3, and a 2 = «3 = 5. The analogue of j7|) is 

\L | = \U'\ = 4r — a±r — 2a 2 r — 3a^r = (4 — ot\ — 2a 2 — 3a^)r. 

The estimate t{n,3 r ) = 6((5 1 / 4 . 2) n ) now yields t(4r,3 r ) = G((5 1 / 4 • 2) 4r ) = G(5 r • 2 4r ). The 
resulting lower bound is 



N = Q* 

= o* 



where 



g . 2 4 2^ 4 ^ ao,ai,a2 ' as ^ ■ 3 ai ■ 5 a2 • 5 Q3 • 2^ 4— a i~ 2a 2— 303) 

in/4 

g _ 2 8_a i _2a 2-3a3+-H"4(ao>c«i,Cf2,a3) . g<*i . 5O2 . 503 



g _ 28— ai— 2a2— 3a3+-H"4(oo, 01,02,03) . 301 _ ^02 _ ^03 



2r 



o* (a n ) 



1/4 



The values a = 0.23, «i = 0.34, a 2 = 0.29, a 3 = 0.14, yield a = 8.6504. . ., and N = O*(8.6504") = 
0(8.65"). The proof of Theorem[T]is now complete. □ 

Remark. The next step in this approach, using .D(n,4 r ), where n = lOr + 2, seems to bring 
no further improvement in the bound. Omitting the details, it gives (with the obvious extended 
notation) : 



N > O* ( 6 • 2a2 ~ 303—404+^5(00,01,02,03,04) _ ^Qi _ g02 _ ]^4 a 3 . 14 Q 4 



n/5 



O* (a r 



where 



g _ 2 10 - a l~ 2Q 2-3o3-4o4+_ff5(oo,Oi,02,Q3,0 4 ) _ ^Oi _ g0 2 _ ]^4 Q 3 . 14° 4 



1/5 



Numerical experiments suggest that the maximum value of a is smaller than 8.65, with the corre- 
sponding a-tuple near ao = 0.127, ot\ = 0.254, a 2 = 0.286, 0:3 = 0.222, 04 = 0.111, which only 
yields a = 8.6485 .... Similarly, when taking reflex chains one unit longer (using L>(n,4 r ), where 
n = 12r + 2), numerical experiments suggest that the maximum value of a is smaller than 8.64. 

3 Lower bound on the maximum number of non-crossing spanning 
trees and forests 

In this section we derive a new lower bound for the number of non-crossing spanning trees on 
the double-chain Z?(n,O r ), hence also for the maximum number of non-crossing spanning trees an 
n-element planar point set can have. The previous best bound, O(10. 42 n ), is due to Dumitrescu [9]. 
By refining the analysis of [9] we obtain a new bound 0(12. 00 n ). 



2 These values have been obtained by numerical experiments. 
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Theorem 2. For the double chain D(n,O r ), we have 



O(12.00 n ) < st(D(n,O r )) < 0(24.68 n ), and 
0(12.26") < cf(D(n,O r )) < 0(24.68"). 

Consequently, st(ra) = 0(12.00") and cf (n) = 0(12. 26 n ). 

Proof. Similarly to [9J, instead of spanning trees, we count (spanning) forests formed by two trees. 
One of the trees is associated with the lower chain L and is called lower tree, while the other tree 
is associated with the upper chain U and is called upper tree. Since the two trees can be connected 
in at most 0(n 2 ) ways, it suffices to bound from below the number of forests of two such trees. 




Figure 3: A double chain with lower and upper tree and four bridges. 

Fig. [3] shows an example. We count only special kinds of forests: no edge of the lower tree 
connects two vertices of the upper chain, and similarly, no edge of the upper tree connects two 
vertices of the lower chain. We call the connected components of the edges between U and L 
bridges. For the class of forests we consider, bridges are subtrees of the lower or the upper tree. A 
bridge is called an (i, j)-bridge if it has % vertices in L and j vertices in U. Every bridge is part 
of either the upper or the lower tree. In the first case we say that the bridge is oriented upwards 
and in the latter case that it is oriented downwards. Since edges cannot cross, the bridges have a 
natural left-to-right order. Fig. [3] shows four bridges, where the first bridge is an upward oriented 
(2, 2)-bridge. We consider only bridges (i,j), with 1 < i,j < z, for some fixed positive integer z. 
For z = 1, our analysis coincides with the one in [9], and we first re-derive the lower bound of 
17(10.42") given there. Successive improvements are achieved by considering lager values of z. 

Let m = re/2 be the number of points in each chain. The distribution of bridges is specified 
by a set of parameters ctij, to be determined later, where the number of (i, j')-bridges is a^rre. To 
simplify further expressions we introduce the following wildcard notation: 



o, 



^"ifc, a *j = , ^2 a kj, and a** = ^a*fc = ^afc*. 



k=l 



k=l 



k=l 



k=l 



A vertex is called a bridge vertex, if it is part of some bridge, and it is a tree vertex otherwise. We 
denote by a^rre the number of bridge vertices along L, and by ajjm the number of bridge vertices 
along U, we have 



oil = kctk*, and a\j = ^ ka*k- 



fc=i 



fc=i 
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To count the forests we proceed as follows. We first count the distributions of the vertices that 
belong to bridges on the lower (Nl) and upper chain (Njj). We then count the number of different 
ways in which the bridges can be realized (Abridges) an d the number of ways in which the bridges 
can be connected to the two trees (-/Vn n k s ). Finally, we estimate the number of the trees within the 
two chains (Agrees)- All these numbers are parameterized by the variables a,j. 

Consider the feasible locations of bridge vertices at the lower chain. We have (q,"^) ways to 
choose the bridge vertices that are in L. Every bridge vertex belongs to some (i, j)-bridge. The 
vertices of the bridges cannot interleave, thus we can describe the configuration of bridges by a 
sequence of tuples that denotes the appearance of the a**m bridges from left to right on L. 
There are ( w m c ^*^ 1 azzm ) such sequences. This give us a total of 

Nl . = f m \( a ** m ^) = e* (2 H{aL)m+a ** H (z 2 ) iail/a ** , ---' azz/a " )m ) 



a.\\m, ai2m, . . . , a zz m 



such "configurations" of bridge vertices along L. 

We now determine how many options we have to place the bridge vertices on U. Since we have 
already specified the sequence of the (i, j)-bridges at the lower chain, all we can do is to select the 
bridge vertices in U. This gives 



N v :-- 



( m ) = e* ( 2 H ^ m ) 

\aumj \ ) 



possibilities for the configuration on U. 

We now study in how many ways the bridges can be added to the two trees. Since all bridges are 
subtrees, we can link one of the bridge vertices with the lower or upper tree. From this perspective 
the whole bridge acts like a super-node in one of the trees. The orientation of each bridge determines 
the tree it belongs to: upwards bridges to the upper tree, and downwards bridges to the lower tree. 
For every pair we orient half of the (i, j)-bridges upwards and half of them downwards. To 
glue the bridges to the trees we have to specify a vertex that will be linked to one of the trees. 
Depending on the orientation of the (i, j)-bridge, we have i candidates for a downwards oriented 
bridge and j candidates for an upward oriented bridge. In total we have 

Minks := n ( ai ' jm /n) ('" 2 r t = n e * ( 2a " jm ) (ur w2 = ©* (2 Q - m ) n 

i,j 1,3 M 

ways to link the bridges with the trees. 

Until now we have specified which vertices belong to which type of bridges, the orientation of 
the bridges, and the vertex where the bridge will be linked to its tree. It remains to count the 
number of ways to actually "draw" the bridges. Let us consider an (i, j)-bridge. All of its edges 
have to go from L to U and the bridge has to be a tree. The number of such trees equals the 
number of triangulations of a polygon with point set {(k, 0) | < k < i} U {(k, 1) | < k < j}. By 
deleting the edges along the horizontal lines y = and y = 1, we define a bijection between these 
triangulations and the combinatorial types of (i, j)-bridges. The number of triangulations is now 
easy to express similarly to Equation ([5]): We have i + j — 2 triangles, and each triangle is adjacent 
to a horizontal edge along either y = or y = 1, where exactly i — 1 triangles are adjacent to line 
y = 0. In total we have Bij := i^^J^) different triangulations and therefore we can express the 
number of different bridges by 

Abridges = ( IT ^ij J 



in 
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Figure 4: All £34 = 10 combinatorial types of (3, 4)-bridges. If an edge differs from its predecessor at the 
top we write a 0, otherwise a 1. We obtain a bijection between the bridges and sequences with three Is and 
two 0s. 



Observe that the upper and the lower trees are trees on a convex point set. By considering the 
bridges as super-nodes, we treat the lower chain as a convex chain of ni vertices. Similarly, we 
think of the upper chain as a convex chain with nu vertices. We have 

nu = (l~y] — - — a k ^m, and n L = (l - ^ — - — a* fc )m. 

k=l k=l 

(Notice that the bridges take away all of its vertices, except one, depending on the orienta- 
tion.) Since the number of non-crossing spanning trees on an n-element convex point set equals 
0*((27/4)") p3], the number of spanning trees within the two chains is given by 

A r trees = 0*((27/4) ni+ " C7 ). 

To finish our analysis we have to find the optimal parameters aij such that 

st(D(n, r )) = O* (N L ■ Nu ■ iV brid gcs ■ blinks • AWs) (8) 

is maximized. 

Assume that we picked all values, except one a uv , and let g U v(~) be the bound of ((8]) in terms 
of this only unspecified parameter. It can be observed that g uv = g vu . Hence we can replace all 
occurrences of a uv in (|8|) with u > v by a vu . To solve the optimization problem we restrict ourselves 
to a few small values of z. Let us start with the case z = 1. Here our analysis coincides with the one 
in [9]. We have only one parameter an, which maximizes ([8]) at 4/(4 + 3\/6) = 0.35 . . . and yields 
a lower bound of 17(10.42™). A first improvement is achieved by considering z = 2. In this case ([8]) 
is bounded by 0(11.61"'), which is attained at an = 0.18, a\i = 0.055, 022 = 0.032. For the next 
step z = 3, we obtain a lower bound of 0(11.89"'), attained at an = 0.15, ai2 = 0.043, 013 = 0.010, 
022 = 0.023, a23 = 0.0085, and 033 = 0.0040. The optimal solution for z = 4 are listed in Table [3j 
The induced bound equals 0(12.00"). By further increasing z we get improved bounds. The best 
bound was computed for z = 8, namely 0(12.0026"). The parameters realizing this bound can be 
found in Table [H For larger z we were not able to complete the numeric computations and got 
stuck at a local maximum, whose induced bound was smaller than 0(12.0026"). All optimization 
problems were solved with computer algebra software. 

Cycle-free graphs. The same approach can be used to bound the number of cycle-free graphs 
(i. e., forests) on the double chain. To this end, we update Equation ([8]) by substituting the quantity 

Agrees with 

^forests = ^* (8.22469"^+"^) , 

which is obtained from the bound for the number of forests on a convex point set p3]. Since we 
picked the bridges such that they introduce no cycles, all graphs considered by our scheme are 
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i = i 


j = 2 


J = 3 


j = A 




3=± 3=2 


J = 3 


j = 4 


i = l 


0.149 


0.0403 


0.00945 


0.00208 


i = l 


0.151 0.0382 


0.00835 


0.00172 


i = 2 




0.0218 


0.00767 


0.00226 


i = 2 


0.0192 


0.00632 


0.00226 


i = 3 






0.00359 


0.00132 


i = 3 




0.00276 


0.00094 


i=A 








0.00058 


i = A 






0.00039 



Table 3: On the left: Parameters that attain the maximum 0(12.00") for st(_D(n,0 r )) for z = 4. On the 
right: Parameters that attain the maximum 0(12.26") for cf (D(n, r )) for z = 4. 



a^ 


i = 1 i=2 


i = 3 


i = 4 


i = 5 


i = 6 


i = 7 


i = 8 






3 = 


1 


0.144 0.0389 


0.00908 


0.001994 


0.000422 


0.0000856 


0.0000152 


1.76 


10" 


-6 


3 = 


2 


0.0209 


0.00733 


0.00214 


0.000569 


0.000140 


0.0000313 


5.12 


10" 


-6 


3 = 


3 




0.00342 


0.00125 


0.000397 


0.000113 


0.0000290 


5.50 


10" 


-6 


3 = 


4 






0.000548 


0.000202 


0.0000655 


0.0000181 


3.34 


10" 


-6 


3 = 


5 








0.0000845 


0.0000298 


8.33 • 10~ 6 


1.25 


10" 


-6 


3 = 


6 










0.0000107 


2.44 • 10~ 6 


5.10 


10" 


-7 


3 = 


7 












1.09 • 10~ 6 


1.31 


10" 


-7 


3 = 


8 














6.97 


10- 


-8 



Table 4: Parameters that realize the maximum 0(12.0026") for st(D(n, r )) for z = 8. 



cycle-free graphs. When the super-node (that represents a bridge) is a singleton in the forest of 
the chain, the combined forest could be constructed by different gluings of that bridge. For this 
reason we have to ignore the possibilities that were covered by -/Vvmks to avoid over-counting. Thus 
we end up with 

cf (D(n, r )) = n* {N L ■ N v ■ N hridges ■ AWsts) • 

For z = 2, our method yields the bound 17(11.94") attained at an = 0.18, a\2 = 0.039, and 
a 2 2 = 0.021; and for z = 3 the bound 0(12.16") attained at a n = 0.12, a 12 = 0.031, a 13 = 0.0080, 
ct22 = 0.016, «23 = 0.0061, and 033 = 0.0031. The bound for z = 4, namely 0(12.26"), is attained 
at the parameters listed in Table El The best bound was obtained for z = 9. In this case the 
parameters in Table [5] imply a bound of 0(12.2618"). Numeric computations for larger numbers 
of z gave no improvement on the bound for cf (D(n, r )). The previous best lower bound by 
Aichholzer et al. was 0(11.62"). 



a^ 


i = 1 


i=2 i = 3 


i = 4 


i = 5 


i = 6 


i = 7 


i = 8 


i = 


9 




3 = 


1 


0.11 


0.028 0.0069 


0.0017 


0.00042 


0.00010 


0.000024 


5.3 • 10^ 


1.2 


• 10" 


-6 


3 = 


2 




0.014 0.0051 


0.0017 


0.00052 


0.00015 


0.000042 


0.000011 


2.7 


• 10- 


-6 


3 = 


3 




0.0025 


0.0010 


0.00038 


0.00013 


0.000042 


0.000012 


3.4 


• 10- 


-6 


3 = 


4 






0.00051 


0.00022 


0.000086 


0.000030 


0.000010 


3.1 


• 10" 


-6 


3 = 


5 








0.00011 


0.000047 


0.000018 


6.5 • 10~ 6 


2.2 


• 10" 


-6 


3 = 


6 










0.000024 


9.4- 10~ 6 


3.7- 10~ 6 


1.4 


• 10" 


-6 


3 = 


7 












5.6- 10~ 6 


1.8 • 10~ 6 


7.5 


• 10" 
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1.5 • 10~ 6 
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• 10" 


-7 
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4.4 
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0(12.0026") for cf (D(n 
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z = 1 


z = 2 


z = 3 


z = 4 


z = 5 


z = 6 


z = 7 


z = 8 


z = 9 


St 


10.424 


11.611 


11.899 


12.004 


11.952 


11.998 


12.002 


12.002 




cf 


11.092 


11.944 


12.169 


12.260 


12.251 


12.258 


12.260 


12.261 


12.261 



Table 6: Bases of the asymptotic exponential bounds for the number of spanning trees st and the number 
of cycle-free graphs cf for the point set D(n, 1 *)) in terms of the parameter z. 

Notice that our analysis relies on numeric computations. We were able to compute the optimal 
parameters up to z = 8 (for spanning trees) and up to z = 9 (for cycle-free graphs). Table [6] lists 
the computed bounds with respect to the chosen parameter z. No significant improvements can be 
expected by solving the induced maximizations problems for larger values of z. 

For comparison, we compute an easy upper bound for both st(D(n, r )) and cf (D(n, r )). 
Every forest F on the double chain splits into three parts: the induced subgraph on L {F£), the 
induced subgraph on U (Fjj) and the remaining part given by the bridge edges (Fb). The graphs 
Fl, Fu, and Fb are non-crossing forests. Since the number of forests on n points in convex position 
is 0(8.225"), there are 0(8.225 n ) different graphs F^UFjj (we rounded up to avoid the O* notation). 
It remains to bound the number of graphs Fb- Every graph counted in Fb can be turned into a 
triangulation by deleting the "non-bridge vertices" and adding the appropriate edges on L and U. 
Let k be the number of bridge vertices on each chain. We have (Tj ways to select the bridge 
vertices. By Equation ([5]) we have 0*(4 fc ) triangulations of the bridge vertices. In total we can 
bound the number of graphs Fb from above by Y^T=q (T) ^ ■ ^° bound the exponential growth 
of this sum we find the dominating summand. Notice that ()I/) 2 4 fc ' is maximized when (^/)2 fc ' is 
maximized. Since ££L (™)2 fc = 3 m , we have (™)2 k ' < 3 m , and so 

(Z) V <*" = *•■ 

Therefore the number of cycle-free non-crossing graphs on D(n, r ) is bounded from above by 
0(ra-3 n -8.225 n ) = 0(24. 68 n ), which is also an upper bound for the number of non-crossing spanning 
trees on D(n, r ). □ 

4 Upper bound for the number of non-crossing spanning cycles 

Newborn and Moser [20J asked what is the maximum number of non-crossing spanning cycles for n 
points in the plane, and they proved that fi((10 1//3 ) n ) < sc(n) < 0(6 n |_§ J!)- The first exponential 
upper bound sc(n) < 10 13n , obtained by Ajtai et al. [2], has been followed by a series of improved 
bounds, e. g., see [6l [U [24]; a more comprehensive history can be found in [7]. The current best 
lower bound sc(ra) > 4.462 n is due to Garcia et al. [15]. The previous best upper bound O(70. 22 n ) 
follows from combining the upper bound 30 n / 4 of Buchin et al. [6j on the number of spanning 
cycles in a triangulation with a new upper bound of tr(n) < 30 n on the number of triangulations 
by Sharir and Sheffer |26] (30 5 / 4 = 70.21 . . .). The bound by Buchin et al. [6] cannot be improved 
much further, since they also present triangulations with 0(2.0845") spanning cycles. However, 
the approach of multiplying tr(n) with the maximum number of spanning cycles in a triangulation 
seems rather weak, since it potentially counts some spanning cycles many times. 

Let C be a spanning cycle. We say that C has a support of x, and write supp(G) = x if C is 
contained in x triangulations of the point set. 
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Observe that a spanning cycle C will be counted supp(C) times in the preceding bound. To 
overcome this inefficiency, we use the concept of pseudo-simultaneously flippable edges (ps-flippable 
edges, for short), introduced in [16] . A set F of edges in a triangulation is ps-flippable if after 
deleting all edges in F, the bounded faces are convex and jointly tile the convex hull of the points. 
One can obtain a lower bound for the support of a spanning cycle C in terms of the number of 
ps-flippable edges that are not in C. The following two properties of ps-flippable edges derived 
in [16] are relevant to us: 



(i) Every triangulation on n points contains at least § — 2 ps-flippable edges. 

(ii) Let C be a spanning cycle. Consider a triangulation T that contains C, and has a set F of 
ps-flippable edges. If j or more edges of F are not contained in C, then supp(C) > 2 3 . 

We now use ps-flippable edges to improve the upper bound for the number of spanning cycles. 

Theorem 3. sc(n) = O (68.62 n ) . 

Proof. Let P be a set of n points in the plane. Assume first that n is even. The exact value of 
sc(P) is 

where the first sum is taken over all triangulations of P, and the second sum is taken over all 
spanning cycles contained in T. 

Consider a triangulation T over the point set P. By property ([!]), there is a set F of n/2 — 2 
ps-flippable edges in T. We wish to bound the number of spanning cycles that are contained in T 
and use exactly k edges from F . The support of any such cycle is at least 2' F '~ fc by property (ii). 
There are ('^') < ( n ^ 2 ) ways to choose k edges from F. 

Since n is even, every spanning cycle in T is the union of two disjoint perfect matchings in T. 
Instead of spanning cycles containing exactly k edges of F, we count the number of pairs of disjoint 
perfect matchings, M\ and M2, that jointly contain exactly k edges of F . Assume, without loss of 
generality, that M2 contains at least as many edges from F as M\. That is, Mi contains at most 
k/2 edges from F, and M2 at least k/2. 

As observed by Buchin et al. [6], a plane graph with v vertices and e edges contains at most 
0{{2e/v) v ^) perfect matchings. By Euler's formula, a triangulation on n points has at most 3n — 6 
edges. Moreover, we know that both, M\ and M2, have to avoid the \F\ — k edges from F that 
were not selected, thus it is safe to remove \F\ — k edges from the triangulation. Since 

(3n - 6) - (n/2 - 2 - k) = 5n/2 + k - 4 < 5n/2 + k, 

the remaining graph has at most 5n/2 + k edges and by the previous expression of Buchin et al. , 
there are O ((5 + 2fc/n)(™/ 4 )) ways to choose Mi. 

After choosing Mi, we know of at least k/2 edges of F which participate in M2 (the edges that 
were chosen from F and do not participate in M\). We can remove the endpoints of these edges 
(together with all incident edges). The resulting graph has at most n — k vertices. We can also 
remove all remaining edges of F and Mi, since they cannot be in M2. It can be easily checked that 
overall we removed n — 2 — k/2 edges (recall that \M\\ = n/2, and that it can hold at most k/2 
edges of F). Thus, given a matching Mi, the number of candidates for M2 is 

// 2(3n-(n-fc/2)) \ (n - fc)/4 \ = Q (fin + k 



n — k J / \ \ n — k 
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We conclude that the number of spanning cycles containing exactly k edges of F is 



° (T)- (5 + 2A;/n)(n/4) -( : 



n — k 



When k is small, i. e., k < an, for some a 6 (0, 1/2), it is better to use the bound 30 n//4 instead 
(which bounds the total number of spanning cycles in T). Substituting these bounds into Q yields 

(go + | o ((f) ■ (^P* ; 

The maximum is attained for a ~ 0.466908, which implies 

sc(P) = O (2.28728 71 ) = tr(P) • O (2.28728 n ) . 

T 

The upper bound sc(n) = 0(68.62") is immediate by substituting the upper bound tr(P) < 30 n 
from [26] into the above expression. 

Finally, the case where n is odd can be handled as follows. Create a new point p outside of 
conv(P), and put P' = PU {p}. It can be shown that sc(P) < sc(P') by mapping every spanning 
cycle of P to a distinct spanning cycle of P' . Given a non-crossing cycle C of P, p can be connected 
to the two endpoints of some edge of C without crossings [171 Lemma 2.1]. We can then apply the 
above analysis for P' , obtaining the same asymptotic bound. □ 



5 Weighted geometric graphs 

In this section we derive bounds on the maximum multiplicity of various weighted geometric graphs 
(weighted by Euclidean length). 



5.1 Longest perfect matchings 

Let n be even, and consider perfect matchings on a set of n points in the plane. It is easy to 
construct n-element point sets (no three of which are collinear) with an exponential number of 
longest perfect matchings: [9] gives constructions with S7(2 n / 4 ) such matchings. Moreover, the 
same lower bound can be achieved with yet another restriction, convex position, imposed on the 
point set; see [9]. Here, we present constructions with an exponential number of maximum (longest) 
non-crossing perfect matchings. 

Theorem 4. For every even n, there exist n-element point sets with at least 2 L^/^J longest non- 
crossing perfect matchings. Consequently, pm max (n) = J7(2™/ 4 ). 

Proof. Assume first that re is a multiple of 4. Let £4 = {a, b, c, d} be a 4-element point set such that 
segment ab is vertical, col lies on the orthogonal bisector of ab (hence, \ac\ = \bc\ and \ad\ = \bd\), 
\ab\ = \cd\ = — and min{|ac|, \ad\} = \ac\ = \bc\ = 2n. Then S4 has two maximum matchings, 
{ac,bd} and {ad, be}, each of which has length at least An. Let the n-element point set P be the 
union of re/4 translated copies of S4 lying in disjoint horizontal strips such that the copies of a are 
almost collinear, all the copies of points a and b lie in a disk of unit diameter, and all the copies of 
points c and d lie in a disk of unit diameter; see Fig. G3 



15 



Figure 5: Left: Two possible maximum matchings for the point set S4 = {a, 6, c, d}. Right: A set of n = 16 
points that admit 2 4 maximum non-crossing perfect matchings. 

If we combine the maximum matchings of all copies of S4, then we obtain 2 n / 4 non-crossing 
perfect matchings of P. All these matchings have the same length, which is at least j • 4n = n 2 . 
We show that this is the maximum possible length of a non-crossing perfect matching of P. Let L 
be the set of the translated copies of points a and 6; and let R be the set of the translated copies 
of points c and d. Then the diameter of L (resp., R) is at most 1, by construction. The length 
of any edge between L and R is at most 2n + 1, while all other edges have length at most 1. If a 
matching has k edges between L and R, then its length is at most k(2n + 1) + — k) = 2kn + ^, 
which is less than n 2 for k < n/2. Hence a maximum matching on P is a bipartite graph between 
L and R. To avoid crossings, every edge in such a bipartite graph must connect points in the same 
copy of 54. 

If n is a multiple of 4 plus 2, the above construction is modified by adding a pair of points at 
distance 2n, placed horizontally above all copies of S4. This pair of points must be matched in 
any maximum non-crossing perfect matching of P. In both cases, P admits at least 2L n / 4 J longest 
non-crossing perfect matchings, as required. □ 

5.2 Non-crossing spanning trees 

For minimum spanning trees, [9] gives constructions for n-element point sets that admit f2(2 n / 4 ) 
such trees, and moreover, this lower bound can be attained with points in convex position. All 
these constructions give non-crossing spanning trees. For maximum spanning trees, [9] gives the 
following construction: start with two points, a and b, and suitably place the remaining n — 2 points 
on the perpendicular bisector of segment ab. While this configuration admits f2(2 n ) maximum non- 
crossing spanning trees (which have maximum weight over all spanning trees), it uses (a large 
number of) collinear points. Next we show that an exponential bound, 0(2 n ), can be achieved 
without allowing collinear points, and moreover, with points in convex position. 

Theorem 5. The vertex set of a regular convex n-gon admits $7(2™ ) longest non-crossing spanning 
trees. Consequently, st max (n) = S7(2 n ). 

Before proving the theorem we introduce some notation, which we will also use in Section El For 
a point set in convex position the span of an edge xy is the smallest number of convex hull edges 
one has to traverse when going from x to y along the convex hull (in clockwise or counterclockwise 
direction). The weight of a tree T, denoted as L(T), is the sum of its edge lengths. 

Proof of Theorem [5], Let P be the vertex set of a regular n-gon inscribed in a circle of unit 
radius. Let p £ P be an arbitrary element of P, and let S p be the star centered at p (consisting of 
segments connecting p to the other n — 1 points). We first prove a counterpart for spanning trees 
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(Lemma [T]) of Alon et al. [3j Lemma 2.2], established for non-crossing matchings. Our argument is 
similar to that used in [3]. 

Lemma 1. A longest non-crossing spanning tree on P has weight L(S P ). 

Proof. Write n = 2k if n is even, or n = 2k + 1 if n is odd. Consider an arbitrary non-crossing 
spanning tree T of P. Notice that the span of any edge is an integer between 1 and k. For 
i = 1, ... ,k, let Ni denote the number of edges of T whose span is at least %. We need the following 
claim. 

Claim. For every i = 1, . . . , k, we have iVj < n — 2i + 1. 

Proof of Claim: The claim trivially holds for N{ < 1, so we can assume that N > 2. One can 
easily check (since T is non-crossing) that there exist two edges P1P2 and P3P4 in T (where j»2 could 
coincide with ^3), each of span at least i, with the following properties: (a) Pi,P2,P3,P4 appear in 
this (clockwise) order on the circle, and (b) no other edge of T with span at least i has an endpoint 
among the points between p\ and P2 or between p% and p^. There are at least i — 1 points on each 
of the two open circular arcs defined by P1P2 and p^Pa- Hence all edges of T with span at least i 
are induced by at most n — 2(i — 1) points. Since this set of edges forms a forest, it has no more 
than n — 2{i — 1) — 1 = n — 2i + 1 elements, as required. □ 

To finalize the proof of the lemma, we show that L(T) < L(S P ). For i = 1, . . . , k, let £i denote 
the (Euclidean) length of an edge with span i. Note that l\ < £2 < ■ ■ ■ < £ki and define Iq = 0. 
Put also iVfc +1 = 0, and observe that the number of edges of span i is equal to N — N{ + \. Consider 
first the case n = 2k + 1. Straightforward algebraic manipulation gives 

k k 

L(T) = J2(Ni - N i+1 )£i = Ni(£i - 4-i) 

i=l i=l 
k k 

< ^2(n -2i + l)(£i - £i_ x ) = 2j2 £ i= L{S p ). 

i=l i=l 

Similarly, for n = 2k we get 

fc-i 

L(T)<2j2ti + Zk = L(S p ). 

i=i 

This concludes the proof of the lemma. □ 

To finalize the proof of Theorem [5j it remains to show that P admits Q,(2 n ) non-crossing 
spanning trees of weight L(S P ). Denote the points in P by pi,P2, ■ ■ ■ ,Pn in clockwise order. We 
construct a family of non-crossing spanning trees on P. Start each tree with the same edge T := 
PiP2- Repeatedly augment T with one edge at a time, such that the clockwise arc spanned by the 
latest edge increases by one, either from its "left" (encode this choice by 0), or from its "right" 
endpoint (encode this choice by 1). Notice that regardless of the choice, the latest edge has always 
the same length. As a result, the clockwise arc spanning the new tree T now includes one more 
point: either its right endpoint moves clockwise by one position, or its left endpoint moves counter- 
clockwise by one position. After n — 2 steps, T is a non-crossing spanning tree of P. If n = 2k + 1, 
the sequence of edge lengths (including the first edge, P1P2) is 

hi ■ ■ ■ > 4-1)4)4-1) • • • ,4- 
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If n = 2k, the sequence of edge lengths (including the first edge, P1P2) is 

^1 > • • • j ^k— 1) *fe) ^k—li ■ ■ ■ > *1 ■ 

In both cases, the weight of each tree is L(S P ). All 0-1 sequences of length n — 2 yield different 
trees, so there are at least 2 n ~ 2 longest non-crossing spanning trees, as required. This completes 
the proof of Theorem [5j □ 



5.3 Longest non-crossing tours 

We show here that the maximum number of longest non-crossing spanning cycles is also exponential 
inn. In the next section we also show that the maximum number of shortest non-crossing spanning 
cycles on n points is exponential in n (Theorem [7]) . 

Theorem 6. Let sc max (n) denote the maximum number of longest non-crossing spanning cycles 
that an n- element point set can have. Then we have Sl(2 n / 3 ) < sc max (n) < 0(68. 62 n ). 

Proof. Since we are counting non-crossing tours, we have sc max (n) = 0(68. 62 n ) by Theorem [3j It 
remains to show the lower bound. For every k G N, we construct a set Q of 4/c + 1 points that 
admits 2 k = 0(2 n / 4 ) longest non-crossing tours. We start by constructing an auxiliary set P of 
2k points. The auxiliary point set P may contain collinear triples, however our final set Q does 
not. Recall that two segments cross if and only if their relative interiors intersect. We construct 
P = {cj, X{ : % = 1, 2, . . . , k} with the following properties: 

(i) for every Xi, the farthest point in P is Cj; 

(ii) the perfect matching M = {ciXi : i = 1, 2, . . . , k} is non-crossing; and 

(iii) the convex hull of P is conv(P) = (x\, c\, c 2 , . . . , c^). 

Note that property (i) implies that M is the maximum matching of P. 

For k G N, let a — T^r. We construct P — {cj,Xj : i — 1,2,..., A;} iteratively. During the 
iterative process, we maintain the properties that 

IxjCjl > max \xicA and |xj + iCj| > max |xj + ic 7 -|. 

Initially, let ci = (0,0), x x = (2,0), and x 2 = (2-^,0). Let l x be a ray emitted by x\ and incident 
to c\. Refer to Fig. EJ If q, Xi and are already defined, we construct points Cj + i and x i+2 (in 
the last iteration, only q+i) as follows. Let £{ + \ be a ray emitted by Xj+i such that Z(£j + i,£j) = a. 
Compute the intersections of ray £{ + i with the circle centered at X{ of radius |xjCj| and the circle 
centered at xi + \ of radius Let q+i G 4+i be the midpoint of the segment between these 

two intersection points. This choice guarantees that |xj + iCj + i| > Ixj+iCjl and |xjCj + i| < \xjCj\ for 
all j < i. Now let G Q+iXj+i be a point at distance at most \ from Xi + \ such that we have 
|xj + 2Cj + i| > |xj + 2Cj| for all j < i. This completes the description of P. 

Note that |xjXj + i| < t- for i = 1, . . . , k — 1, and so the points x±, . . . , Xk lie in a disk of diameter 
1. Hence, for every point Xi, the farthest point in P is in {cj : j = 1, ...,k}. By the above 
construction, the farthest point from Xi in {cj : j = 1, . . . , k} is Cj. This proves that P has property 
(i). It is easy to verify that P has properties (ii) and (iii), as well. 

We now construct the point set Q based on P. Let 5 > be a sufficiently small constant. For 
every segment CjXj we construct a skinny deltoid Aj = (a^, 6j, q, dj), see Fig. [6j such that G CjXj 
is at distance 5 from Xj, we have |6jCj| = |cjdj| = 5, and |aj6j| = |ajCj| = |ajdi| = |cjXj| — 5. 
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Figure 6: Left: The auxiliary point set P for k = 3. Right: A long and skinny deltoid Aj = (dj, bi, Ci, di). 

Since the segments CjXj are pairwise non-crossing and 5 > is small, the deltoids A, are pairwise 
interior disjoint. Let Q be the set of vertices of all deltoids Aj, i = 1, . . . ,k, and the point x\. 
Since conv(P) = (ci,c 2 , . . . , Ck,xi), we have conv(Q) = (&i, ci, d%, 62, C2,d 2 , ■ ■ ■ ,h, Ck,dk,x\), and 
the points {aj : i = 1, . . . , k} lie in the interior of conv(<5). If 5 > is sufficiently small, then the 
farthest points from aj in Q are 6j, Cj, and di, for every i = 1,2,... , A;. 

Every non-crossing tour of Q visits the convex hull vertices in the cyclic order determined by 
conv((5). We obtain a non-crossing tour by replacing some edges of coiw(Q) with non-crossing 
paths visiting the points lying in the interior of conv(Q). If we replace either edge 6jCj or Cjdj with 
the path (pi, aj, a) or (a, aj, di), respectively, for every i = 1, 2, . . . , k, then we obtain a tour. Let 
H be the set of 2 k tours obtained in this way. These tours are non-crossing, since for every i, we 
exchange an edge of Aj with a path lying in Aj, and the deltoids Aj are interior disjoint. The 
tours in % have the same length, L = |conv(Q)| — k5 + 2^^ =1 |ajCj|, since |aj&j| = |ajdj| = |ajCj|. It 
remains to show that this length is maximal. Note that a non-crossing tour cannot have an edge 
between two non-consecutive vertices of conv(Q). Hence, every edge that intersects the interior of 
conv((5) must be incident to some point aj lying in the interior of conv(Q). Each aj is incident to 
two edges of a tour: The total length of these two edges is at most 2|ajCj|, which is attained if aj 
is connected to bi, Ci, or dj. In any other case, it is less than 2|ajCj| — k5 if 5 > is sufficiently 
small. The total length the edges on the boundary of the convex hull is less than |conv(Q)|. So 
any non-crossing tour in which some point a, is not connected to 6j, Cj or to Cj, dj must have length 
less than L. This implies sc max (n) = f2(2 n / 4 ). 

To obtain the asserted bound, we use a skinny hexagon (instead of deltoid Aj) with five equidis- 
tant vertices on a circle centered at aj. We now have four possible ways to insert each aj into the 
tour, which implies sc max (n) = ft(4™/ 6 ) = ft(2"/ 3 ). □ 

Typically for the longest matching, spanning tree or spanning cycle, one expects to see many 
crossings. Somewhat surprisingly, we show that this is not always the case. 

Corollary 1. For every even n>2, there exists an n- element point set (in general position) whose 
longest perfect matching is non- crossing. 

Proof. For every k > 1, consider the set {aj,Cj : i = 1,2,..., A;} of 2k points, and its perfect 
matching {ajCj : i = 1, 2, . . . , k} in the proof of Theorem El This is the longest perfect matching 
(over all perfect matchings of the point set, crossing or non-crossing), and moreover, it is non- 
crossing. □ 
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6 Tours 



In this section we derive estimates on the maximum multiplicity of the shortest (minimum) and, 
respectively, the longest (maximum) Hamiltonian tour on n points (crossings allowed). While the 
shortest tour has the non-crossing attribute for free, the longest tour typically has crossings. 

6.1 Shortest tours 

Theorem 7. Let sc m i n (n) denote the maximum number of shortest tours that an n-element point 
set can have. 

(i) For points in convex position, there is exactly one shortest spanning cycle, i. e., sc m ; n (n) = 1. 

(ii) For points in general position, sc m ; n (n) is exponential in n . More precisely, 2 L^/sj < SCmin ( ra ) < 
0(68.62 n ). 

Proof, (i) It is well known (and easy to prove) that a shortest tour of a convex point set does not 
admit any crossings, thus the vertices have to be visited in clockwise or counterclockwise order. 

(ii) Consider now a set S of n points in general position. Since any minimum tour of S is a 
non-crossing spanning cycle of S, by Theorem [3] we have sc m ; n (n) < sc(n) = 0(68. 62 n ). 

A lower bound of 2 L"-/ 3 J i s given by the following construction. Assume first that n is a multiple 
of 3. Consider an isosceles triangle Aabc with sides e, e and e/4, see Fig. [7] (left). Let S consist of 




Figure 7: Lower bound for shortest tours. 

n/3 rotated copies of A along a circle of unit radius, as in Fig. [7] (right). If e is small enough, the 
shortest tour of S must visit the three vertices in each copy of A in their order along the circle. 
We can assume w.l.o.g. that the groups (copies) are visited in clockwise order. Then the first point 
visited in each group is (a rotated copy of) a. Since each group of three points can be minimally 
traversed in two ways, a, b, c, or a, c, b, the number of shortest tours is 2 n / 3 . The construction can 
be easily modified for any n, by removing one or two points from one of the groups. □ 

6.2 Longest tours for points in convex position 

In the following we give tight bounds for the number of maximum tours on n points in convex 
position (allowing crossings) and outline how to compute such tours efficiently. The result for odd 
n seems to have been first discovered by Ito et al. [19]. Nevertheless we present our own simpler 
proof here for completeness. 
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Theorem 8. Let tc max (n) denote the maximum number of longest tours that an n-element point 
set in convex position can have. For n odd we have tc max (n) = 1 and the (unique) longest tour is 
a thrackle. For n even we have tc max (n) = n/2. 

Proof. Let P be a set of n points in convex position. Assume that n = 2k if n is even, and n = 2k — 1 
if n is odd. We pick an orientation for every possible tour arbitrarily. According to this orientation 
we denote an edge xy of a tour as x — > y if x is the predecessor of y in the tour. We call two edges 
x — > y and u — > v parallel if they are disjoint (including their endpoints) and the segments xu and 
yv are disjoint. For instance, the edges 03 — > a 2 and b 2 — > 63 in Fig. Wb) are parallel. Two disjoint 
edges that are not parallel are called anti-parallel. We say two vertices of P are consecutive if they 
define a convex hull edge. To prove the theorem we first prove three simple observations. 

Observation (a): The longest tour does not contain a pair of anti-parallel edges. 
Suppose the longest tour does contain anti-parallel edges x — > y and u — > v. We delete both edges 
and construct a new tour by adding x — > u and y — > v and re-orienting the part of the old tour 
between u and y. The new tour is longer since in a convex quadrilateral the sum of the diagonals 
exceeds the sum of two opposing side lengths. 

Observation (b): In a longest tour with parallel edges x — > y and u — > v, the vertices x,u and the 
vertices y, v are consecutive. In particular, a longest tour contains no triplet of pairwise parallel 
edges. 

Assume that there exists a vertex w "in between" x and u with predecessor w' . Then the edge 
w' —> w would determine a pair of anti-parallel edges with either x —> y or u —> v, in contradiction 
to Observation (a). By a similar argument one can show that y and v have to be consecutive. 

Observation (c): The span of every edge in a longest tour is at least k — 1. 

Suppose, to the contrary, that the longest tour contains an edge x — > y with span at most k — 2. 
By the pigeonhole principle there is at least one edge u — > v that does not cross x — > y. Due to 
Observation (b), uv has to be an edge of span at least n — k. Let {L, R} be the partition of P\{x, y} 
induced by x — > y (that is, points left and right relative to x — >■ y). Assume further that \R\ > \L\, 
that is, {u, v} C R. By Observation (b), a longest tour cannot contain three pairwise parallel edges. 
Hence, in the longest tour the points in R \ {u} have to have a successor from L U {x}. But, again 
due to the pigeonhole principle, this is impossible, since \R\ — 1 > k — 1, but \L\ + 1 < k — 2. 

We now proceed with the proof of Theorem Assume first that n = 2k — 1 is odd. Due to 
observation (c), the longest tour does not have any edge of span at most k — 2, hence the span 
of every edge is k — 1, which is the largest possible span for n odd. This determines the longest 
geometric tour. For every vertex we have only two possible edges with span k — 1, hence the longest 
tour is unique. 

Now assume that n = 2k is even. Notice that we have to use at least one edge of span k, since 
otherwise (using edges with span k — 1 only) there would be two anti-parallel edges. Assume that a 
longest tour contains the edge a\ — > bi of span k. Let a<i denote the predecessor of a±, and &2 denote 
the successor of b\ in the longest tour. The edges a<i — > a\ and b\ — > b<i need to have span k — 1. 
There are two possibilities for the relative position of 0201 and 6261. In case they cross, all outgoing 
edges from the vertices on the arc between a\ and b\ (in cyclic order of P, opposite from 02 and 
62) have to cross b\ —^62- This yields a contradiction by the pigeonhole principle (see Fig. E^a)). 
Thus we can assume that ai and 62 are located on opposite sides of a\ — > 61, as in Fig. [H^b). We 
extend the partial tour step by step, choosing new edges a^+i — > a^, b{ — > with span k — 1 
appropriately. If we would use an edge with span k we would "close" the spanning cycle before all 
vertices have been visited. We stop after we found n — 1 edges and add the final edge b n / 2 — > o- n /2 
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(a) 



(b) 



Figure 8: (a) A configuration that does not allow an extension to a longest tour, since the possible 5 
"outgoing" edges have only 4 candidates to connect with, (b) A valid extension. The edges a i+1 — > ai and 
b{ — > bi + \ form a pair of parallel edges. 



with span k. To see that this scheme works correctly, notice that the edges added in every step 
form a pair of parallel edges with span k — 1, which "rotates" around P. The construction can be 
characterized by the location of the edges with span k. Since a\ and b n /2 are consecutive, as well 
as a n /2 and b%, we have at most n/2 choices to select these two edges. 

If P is the vertex set of a regular n-gon, then P has n/2 longest tours by rotational symmetry. 
The 5 longest tours for n = 10 are depicted in Fig. [UJ □ 




Figure 9: The 5 longest tours on the vertex set of a regular 10-gon. The two edges with span 5 are drawn 
with bold lines. 

As an immediate algorithmic corollary of Theorem El the longest tours on a set of n points 
in convex position can be computed in 0(n log n) time. If n is odd and the convex polygon 
P = Po,Pi, ■ ■ ■ iPn-i is given, the tour can be computed in 0(n) time: start with i = 0, and 
iteratively set i <— , where the indices are taken modulo n. If n is even, we have n/2 candidates 

for the longest tour. Given a candidate, we can construct the next candidate by exchanging only 4 
edges of the tour. Thus, after computing the weight of the first candidate, all other candidates can 
be evaluated in 0(1) per tour. As in the case of odd n, we can compute the longest tours in 0(n) 
time if the cyclic order of P is given, and in 0(n log n) time otherwise. 

7 Conclusion 

Our investigations leave some questions open: 

1. We used a special point configuration, D(n, 3 r ), to prove a new lower bound on tr(n), 
the maximum number of triangulations. The question arises if this (or a similar) configuration 
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could also give better lower bounds for other quantities we studied, such as cf (n), st(n), or sc(ra). 
However, no analysis for other geometric graph classes has been done so far, not even for the 
simpler-to-analyze double zig-zag chain D(n,l r ),. 

2. While for the maximum number of triangulations, the upper and lower bounds are reasonably 
close, this is not the case for non-crossing spanning cycles. Recall that the current best lower bound, 
offered by the double chain configuration is 17(4.462"), and our new upper bound is only 0(68. 62 n ). 
Most likely, the lower bound is closer to the truth, however the arguments are missing to justify 
this belief. Clearly a large gap remains to be covered. 

3. Is it possible to obtain upper bounds for the maximum multiplicity of weighted configura- 
tions, better than those for the corresponding unweighted structures? For example, is sc max (n) 
significantly smaller than sc(n), or pm min (n) significantly smaller than pm(n)? 

4. For points in convex position, we devised an 0(n log n)-time algorithm for computing the 
longest tour. However, for points in general position it is not known whether this problem is NP- 
hard. Barvinok et al. [4] showed that the problem is solvable in 0(n) time under the L\ metric [4]. 
On the negative side, they also show that the problem is NP-complete under the L2 metric for points 
in 3-space. Nothing is known about the complexity of computing the longest non-crossing spanning 
tour, spanning path, perfect matching or spanning tree. However, constant ratio approximations 
are available for the last three problems; see [3| [12], 

References 

[1] O. Aichholzer, T. Hackl, B. Vogtenhuber, C. Huemer, F. Hurtado and H. Krasser, On the 
number of plane geometric graphs, Graphs and Combinatorics 23(1) (2007), 67-84. 

[2] M. Ajtai, V. Chvatal, M. Newborn and E. Szemeredi, Crossing-free subgraphs, Annals Discrete 
Math. 12 (1982), 9-12. 

[3] N. Alon, S. Rajagopalan, and S. Suri, Long non-crossing configurations in the plane, Funda- 
menta Informaticae 22 (1995), 385-394. A preliminary version in Proc. 9th ACM Symp. on 
Comput. Geom., ACM Press, 1993, pp. 257-263. 

[4] A. I. Barvinok, S. P. Fekete, D. S. Johnson, A. Tamir, G. J. Woeginger, and R. Woodroofe, 
The geometric maximum traveling salesman problem, J. ACM 50(5) (2003), 641-664. 

[5] P. Beame, M. Saks and J. Thathachar, Time-space tradeoffs for branching programs, Proc. 
39th Annual Symposium on Foundations of Computer Science, 1998, pp. 254-263. 

[6] K. Buchin, C. Knauer, K. Kriegel, A. Schulz, and R. Seidel, On the number of cycles in planar 
graphs, Proc. 13th Annual International Conference on Computing and Combinatorics, volume 
4598 of LNCS, Springer, 2007, pp. 97-107. 

[7] E. Demaine, Simple polygonizations, http : / / erikdemaine . org/polygonization/ (version of 
November, 2010). 

[8] A. Dumitrescu, On two lower bound constructions, Proc. 11th Canadian Conf. on Computa- 
tional Geometry, 1999, pp. 111-114. 

[9] A. Dumitrescu, On the maximum multiplicity of some extreme geometric configurations in the 
plane, G combinatorics 12(1) (2002), 5-14. 



23 



[10] A. Dumitrescu and R. Kaye, Matching colored points in the plane: some new results, Compu- 
tational Geometry: Theory and Applications 19(1) (2001), 69-85. 

[11] A. Dumitrescu and W. Steiger, On a matching problem in the plane, Discrete Mathematics 
211(1-3) (2000), 183-195. 

[12] A. Dumitrescu and Cs. D. Toth, Long non-crossing configurations in the plane, Discrete and 
Computational Geometry 44(4) (2010), 727-752. 

[13] A. Dumitrescu, A. Schulz, A. Sheffer and Cs. D. Toth, Bounds on the maximum multiplicity of 
some common geometric graphs, Proc. 28th International Symposium on Theoretical Aspects of 
Computer Science, Leibniz International Proceedings in Informatics (LIPIcs) 9, 2011, pp. 637- 
648. 

[14] P. Flajolet and M. Noy, Analytic combinatorics of non-crossing configurations, Discrete Math- 
ematics 204 (1999), 203-229. 

[15] A. Garcia, M. Noy and A. Tejel, Lower bounds on the number of crossing-free subgraphs of 
Kn i Computational Geometry: Theory and Applications 16(4) (2000), 211-221. 

[16] M. Hoffmann, M. Sharir, A. Sheffer, Cs. D. Toth, and E. Welzl, Counting plane graphs: 
flippability and its applications, Proc. Algorithms and Data Structure Symposium (WADS), 
vol. 6844 of LNCS, Springer, 2011, pp. 524-535. Also at arxiv.org/abs/1012.0591 

[17] F. Hurtado, M. Kano, D. Rappaport, and Cs. D. Toth: Encompassing colored planar straight 
line graphs, Computational Geometry: Theory and Applications 39 (1) (2008), 14-23. 

[18] F. Hurtado and M. Noy, Counting triangulations of almost-convex polygons, Ars Combinatoria 
45 (1997), 169-179. 

[19] H. Ito, H. Uehara, and M. Yokoyama, Lengths of tours and permutations on a vertex set of a 
convex polygon, Discrete Applied Mathematics 115(1-3) (2001), 63-71. 

[20] M. Newborn and W. O. J. Moser, Optimal crossing-free Hamiltonian circuit drawings of K n , 
J. Combin. Theory Ser. B 29 (1980), 13-26. 

[21] J. Pach and P. K. Agarwal, Combinatorial Geometry, John Wiley, New York, 1995. 

[22] A. Razen, J. Snoeyink, and E. Welzl, Number of crossing-free geometric graphs vs. triangula- 
tions, Electronic Notes in Discrete Mathematics 31 (2008), 195-200. 

[23] F. Santos and R. Seidel, A better upper bound on the number of triangulations of a planar 
point set, Journal of Combinatorial Theory Ser. A 102(1) (2003), 186-193. 

[24] M. Sharir and E. Welzl, On the number of crossing-free matchings, cycles, and partitions, 
SIAM J. Comput. 36(3) (2006), 695-720. 

[25] M. Sharir and E. Welzl, Random triangulations of point sets, Proc. 22nd Annual ACM-SIAM 
Symposium on Computational Geometry, ACM Press, 2006, pp. 273-281. 

[26] M. Sharir and A. Sheffer, Counting triangulations of planar point sets, Electronic Journal of 
Combinatorics 18 (2011), P70. 



24 



